A model industrial workhorse: Bacillus subtilis strain 168 and its genome after a quarter of a century.

The vast majority of genomic sequences are automatically annotated using various software programs. The accuracy of these annotations depends heavily on the very few manual annotation efforts that combine verified experimental data with genomic sequences from model organisms. Here, we summarize the updated functional annotation of Bacillus subtilis strain 168, a quarter century after its genome sequence was first made public. Since the last such effort 5 years ago, 1168 genetic functions have been updated, allowing the construction of a new metabolic model of this organism of environmental and industrial interest. The emphasis in this review is on new metabolic insights, the role of metals in metabolism and macromolecule biosynthesis, functions involved in biofilm formation, features controlling cell growth, and finally, protein agents that allow class discrimination, thus allowing maintenance management, and accuracy of all cell processes. New 'genomic objects' and an extensive updated literature review have been included for the sequence, now available at the International Nucleotide Sequence Database Collaboration (INSDC: AccNum AL009126.4).

A novel binary pesticidal protein from Chryseobacterium arthrosphaerae controls western corn rootworm by a different mode of action to existing commercial pesticidal proteins.

The western corn rootworm (WCR) Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae) remains one of the economically most important pests of maize (Zea mays) due to its adaptive capabilities to pest management options. This includes the ability to develop resistance to some of the commercial pesticidal proteins originating from different strains of Bacillus thuringiensis. Although urgently needed, the discovery of new, environmentally safe agents with new modes of action is a challenge. In this study we report the discovery of a new family of binary pesticidal proteins isolated from several Chryseobacterium species. These novel binary proteins, referred to as GDI0005A and GDI0006A, produced as recombinant proteins, prevent growth and increase mortality of WCR larvae, as does the bacteria. These effects were found both in susceptible and resistant WCR colonies to Cry3Bb1 and Cry34Ab1/Cry35Ab1 (reassigned Gpp34Ab1/Tpp35Ab1). This suggests GDI0005A and GDI0006A may not share the same binding sites as those commercially deployed proteins and thereby possess a new mode of action. This paves the way towards the development of novel biological or biotechnological management solutions urgently needed against rootworms.

Analyzing Prokaryotic Transcriptomics in the Light of Genome Data with the MicroScope Platform.

Large-scale genome sequencing and the increasingly massive use of high-throughput approaches produce a vast amount of new information that completely transforms our understanding of thousands of microbial species occurring in our environment. However, despite the development of powerful bioinformatics approaches, full interpretation of the content of these genomes remains a difficult task. To address this challenge, the MicroScope platform has been developed. It is an integrated Web platform for management, annotation, comparative analysis, and visualization of microbial genomes ( https://mage.genoscope.cns.fr/microscope ). Launched in 2005, the platform has been under continuous development and provides analyzes for complete and ongoing genome projects together with metabolic network reconstruction and transcriptomic experiments allowing users to improve the understanding of gene functions. MicroScope platform is widely used by microbiologists from academia and industry all around the world for collaborative studies and expert annotation. It enables collaborative work in a rich comparative genomic context and improves community-based curation efforts. Here, we describe the protocol to follow for the integration and analysis of transcriptomics data in the Microscope platform. The chapter reviews each key step from the experimental design to the analysis and interpretation of the experiment data and results. The integration of transcriptomics data gives a dynamic view of the genome by allowing the users to improve the understanding of gene functions by interpreting them in the light of regulatory cell processes. Moreover, they can also contribute to the refinement of genome annotation through the discovery of new genes and help to fill metabolic gaps.

Complete Genome Sequence of Tepidibacter sp. Strain 8C15b, Isolated from Bank Sediments of Haiphong Bay, Vietnam.

We report the complete genome sequence of Tepidibacter sp. strain 8C15b, isolated from bank sediments of Haiphong Bay, Vietnam. The genome includes a 3,628,320-bp circular chromosome and a plasmid of 38,213 bp.

Correction: PPanGGOLiN: Depicting microbial diversity via a partitioned pangenome graph.

[This corrects the article DOI: 10.1371/journal.pcbi.1007732.].

From Strain Characterization to Field Authorization: Highlights on Bacillus velezensis Strain B25 Beneficial Properties for Plants and Its Activities on Phytopathogenic Fungi.

Agriculture is in need of alternative products to conventional phytopharmaceutical treatments from chemical industry. One solution is the use of natural microorganisms with beneficial properties to ensure crop yields and plant health. In the present study, we focused our analyses on a bacterium referred as strain B25 and belonging to the species Bacillus velezensis (synonym B. amyloliquefaciens subsp. plantarum or B. methylotrophicus), a promising plant growth promoting rhizobacterium (PGPR) and an inhibitor of pathogenic fungi inducing crops diseases. B25 strain activities were investigated. Its genes are well preserved, with their majority being common with other Bacillus spp. strains and responsible for the biosynthesis of secondary metabolites known to be involved in biocontrol and plant growth-promoting activities. No antibiotic resistance genes were found in the B25 strain plasmid. In vitro and in planta tests were conducted to confirm these PGPR and biocontrol properties, showing its efficiency against 13 different pathogenic fungi through antibiosis mechanism. B25 strain also showed good capacities to quickly colonize its environment, to solubilize phosphorus and to produce siderophores and little amounts of auxin-type phytohormones (around 13,051 µg/mL after 32 h). All these findings combined to the fact that B25 demonstrated good properties for industrialization of the production and an environmental-friendly profile, led to its commercialization under market authorization since 2018 in several biostimulant preparations and opened its potential use as a biocontrol agent.

Phylogroup stability contrasts with high within sequence type complex dynamics of Escherichia coli bloodstream infection isolates over a 12-year period.

BACKGROUND: Escherichia coli is the leading cause of bloodstream infections, associated with a significant mortality. Recent genomic analyses revealed that few clonal lineages are involved in bloodstream infections and captured the emergence of some of them. However, data on within sequence type (ST) population genetic structure evolution are rare. METHODS: We compared whole genome sequences of 912 E. coli isolates responsible for bloodstream infections from two multicenter clinical trials that were conducted in the Paris area, France, 12 years apart, in teaching hospitals belonging to the same institution ("Assistance Publique-Hôpitaux de Paris"). We analyzed the strains at different levels of granularity, i.e., the phylogroup, the ST complex (STc), and the within STc clone taking into consideration the evolutionary history, the resistance, and virulence gene content as well as the antigenic diversity of the strains. RESULTS: We found a mix of stability and changes overtime, depending on the level of comparison. Overall, we observed an increase in antibiotic resistance associated to a restricted number of genetic determinants and in strain plasmidic content, whereas phylogroup distribution and virulence gene content remained constant. Focusing on STcs highlighted the pauci-clonality of the populations, with only 11 STcs responsible for more than 73% of the cases, dominated by five STcs (STc73, STc131, STc95, STc69, STc10). However, some STcs underwent dramatic variations, such as the global pandemic STc131, which replaced the previously predominant STc95. Moreover, within STc131, 95 and 69 genomic diversity analysis revealed a highly dynamic pattern, with reshuffling of the population linked to clonal replacement sometimes coupled with independent acquisitions of virulence factors such as the pap gene cluster bearing a papGII allele located on various pathogenicity islands. Additionally, STc10 exhibited huge antigenic diversity evidenced by numerous O:H serotype/fimH allele combinations, whichever the year of isolation. CONCLUSIONS: Altogether, these data suggest that the bloodstream niche is occupied by a wide but specific phylogenetic diversity and that highly specialized extra-intestinal clones undergo frequent turnover at the within ST level. Additional worldwide epidemiological studies overtime are needed in different geographical and ecological contexts to assess how generalizable these data are.

Complete Genome Sequences of Two Pseudomonas Species Isolated from Marine Environments of the Pacific Ocean.

Pseudomonas marincola YsY11 and Pseudomonas oleovorans T9AD were both isolated from marine environments of the Pacific Ocean. Here, we report the whole-genome sequences of these two organisms. Pseudomonas marincola YsY11 consists of a single 4.77-Mb chromosome, and Pseudomonas oleovorans T9AD consists of a 5.57-Mb chromosome and a 2.8-kb plasmid.

panRGP: a pangenome-based method to predict genomic islands and explore their diversity.

MOTIVATION: Horizontal gene transfer (HGT) is a major source of variability in prokaryotic genomes. Regions of genome plasticity (RGPs) are clusters of genes located in highly variable genomic regions. Most of them arise from HGT and correspond to genomic islands (GIs). The study of those regions at the species level has become increasingly difficult with the data deluge of genomes. To date, no methods are available to identify GIs using hundreds of genomes to explore their diversity. RESULTS: We present here the panRGP method that predicts RGPs using pangenome graphs made of all available genomes for a given species. It allows the study of thousands of genomes in order to access the diversity of RGPs and to predict spots of insertions. It gave the best predictions when benchmarked along other GI detection tools against a reference dataset. In addition, we illustrated its use on metagenome assembled genomes by redefining the borders of the leuX tRNA hotspot, a well-studied spot of insertion in Escherichia coli. panRPG is a scalable and reliable tool to predict GIs and spots making it an ideal approach for large comparative studies. AVAILABILITY AND IMPLEMENTATION: The methods presented in the current work are available through the following software: https://github.com/labgem/PPanGGOLiN. Detailed results and scripts to compute the benchmark metrics are available at https://github.com/axbazin/panrgp_supdata.

Construction of a complete set of Neisseria meningitidis mutants and its use for the phenotypic profiling of this human pathogen.

The bacterium Neisseria meningitidis causes life-threatening meningitis and sepsis. Here, we construct a complete collection of defined mutants in protein-coding genes of this organism, identifying all genes that are essential under laboratory conditions. The collection, named NeMeSys 2.0, consists of individual mutants in 1584 non-essential genes. We identify 391 essential genes, which are associated with basic functions such as expression and preservation of genome information, cell membrane structure and function, and metabolism. We use this collection to shed light on the functions of diverse genes, including a gene encoding a member of a previously unrecognised class of histidinol-phosphatases; a set of 20 genes required for type IV pili function; and several conditionally essential genes encoding antitoxins and/or immunity proteins. We expect that NeMeSys 2.0 will facilitate the phenotypic profiling of a major human bacterial pathogen.

Alterocin, an Antibiofilm Protein Secreted by Pseudoalteromonas sp. Strain 3J6.

We sought to identify and study the antibiofilm protein secreted by the marine bacterium Pseudoalteromonas sp. strain 3J6. The latter is active against marine and terrestrial bacteria, including Pseudomonas aeruginosa clinical strains forming different biofilm types. Several amino acid sequences were obtained from the partially purified antibiofilm protein, named alterocin. The Pseudoalteromonas sp. 3J6 genome was sequenced, and a candidate alt gene was identified by comparing the genome-encoded proteins to the sequences from purified alterocin. Expressing the alt gene in another nonactive Pseudoalteromonas sp. strain, 3J3, demonstrated that it is responsible for the antibiofilm activity. Alterocin is a 139-residue protein that includes a predicted 20-residue signal sequence, which would be cleaved off upon export by the general secretion system. No sequence homology was found between alterocin and proteins of known functions. The alt gene is not part of an operon and adjacent genes do not seem related to alterocin production, immunity, or regulation, suggesting that these functions are not fulfilled by devoted proteins. During growth in liquid medium, the alt mRNA level peaked during the stationary phase. A single promoter was experimentally identified, and several inverted repeats could be binding sites for regulators. alt genes were found in about 30% of the Pseudoalteromonas genomes and in only a few instances of other marine bacteria of the Hahella and Paraglaciecola genera. Comparative genomics yielded the hypothesis that alt gene losses occurred within the Pseudoalteromonas genus. Overall, alterocin is a novel kind of antibiofilm protein of ecological and biotechnological interest.IMPORTANCE Biofilms are microbial communities that develop on solid surfaces or interfaces and are detrimental in a number of fields, including for example food industry, aquaculture, and medicine. In the latter, antibiotics are insufficient to clear biofilm infections, leading to chronic infections such as in the case of infection by Pseudomonas aeruginosa of the lungs of cystic fibrosis patients. Antibiofilm molecules are thus urgently needed to be used in conjunction with conventional antibiotics, as well as in other fields of application, especially if they are environmentally friendly molecules. Here, we describe alterocin, a novel antibiofilm protein secreted by a marine bacterium belonging to the Pseudoalteromonas genus, and its gene. Alterocin homologs were found in about 30% of Pseudoalteromonas strains, indicating that this new family of antibiofilm proteins likely plays an important albeit nonessential function in the biology of these bacteria. This study opens up the possibility of a variety of applications.

PPanGGOLiN: Depicting microbial diversity via a partitioned pangenome graph.

The use of comparative genomics for functional, evolutionary, and epidemiological studies requires methods to classify gene families in terms of occurrence in a given species. These methods usually lack multivariate statistical models to infer the partitions and the optimal number of classes and don't account for genome organization. We introduce a graph structure to model pangenomes in which nodes represent gene families and edges represent genomic neighborhood. Our method, named PPanGGOLiN, partitions nodes using an Expectation-Maximization algorithm based on multivariate Bernoulli Mixture Model coupled with a Markov Random Field. This approach takes into account the topology of the graph and the presence/absence of genes in pangenomes to classify gene families into persistent, cloud, and one or several shell partitions. By analyzing the partitioned pangenome graphs of isolate genomes from 439 species and metagenome-assembled genomes from 78 species, we demonstrate that our method is effective in estimating the persistent genome. Interestingly, it shows that the shell genome is a key element to understand genome dynamics, presumably because it reflects how genes present at intermediate frequencies drive adaptation of species, and its proportion in genomes is independent of genome size. The graph-based approach proposed by PPanGGOLiN is useful to depict the overall genomic diversity of thousands of strains in a compact structure and provides an effective basis for very large scale comparative genomics. The software is freely available at https://github.com/labgem/PPanGGOLiN.

MicroScope: an integrated platform for the annotation and exploration of microbial gene functions through genomic, pangenomic and metabolic comparative analysis.

Large-scale genome sequencing and the increasingly massive use of high-throughput approaches produce a vast amount of new information that completely transforms our understanding of thousands of microbial species. However, despite the development of powerful bioinformatics approaches, full interpretation of the content of these genomes remains a difficult task. Launched in 2005, the MicroScope platform (https://www.genoscope.cns.fr/agc/microscope) has been under continuous development and provides analysis for prokaryotic genome projects together with metabolic network reconstruction and post-genomic experiments allowing users to improve the understanding of gene functions. Here we present new improvements of the MicroScope user interface for genome selection, navigation and expert gene annotation. Automatic functional annotation procedures of the platform have also been updated and we added several new tools for the functional annotation of genes and genomic regions. We finally focus on new tools and pipeline developed to perform comparative analyses on hundreds of genomes based on pangenome graphs. To date, MicroScope contains data for >11 800 microbial genomes, part of which are manually curated and maintained by microbiologists (>4500 personal accounts in September 2019). The platform enables collaborative work in a rich comparative genomic context and improves community-based curation efforts.

Ectosymbiotic bacteria at the origin of magnetoreception in a marine protist.

Mutualistic symbioses are often a source of evolutionary innovation and drivers of biological diversification1. Widely distributed in the microbial world, particularly in anoxic settings2,3, they often rely on metabolic exchanges and syntrophy2,4. Here, we report a mutualistic symbiosis observed in marine anoxic sediments between excavate protists (Symbiontida, Euglenozoa)5 and ectosymbiotic Deltaproteobacteria biomineralizing ferrimagnetic nanoparticles. Light and electron microscopy observations as well as genomic data support a multi-layered mutualism based on collective magnetotactic motility with division of labour and interspecies hydrogen-transfer-based syntrophy6. The guided motility of the consortia along the geomagnetic field is allowed by the magnetic moment of the non-motile ectosymbiotic bacteria combined with the protist motor activity, which is a unique example of eukaryotic magnetoreception7 acquired by symbiosis. The nearly complete deltaproteobacterial genome assembled from a single consortium contains a full magnetosome gene set8, but shows signs of reduction, with the probable loss of flagellar genes. Based on the metabolic gene content, the ectosymbiotic bacteria are anaerobic sulfate-reducing chemolithoautotrophs that likely reduce sulfate with hydrogen produced by hydrogenosome-like organelles6 underlying the plasma membrane of the protist. In addition to being necessary hydrogen sinks, ectosymbionts may provide organics to the protist by diffusion and predation, as shown by magnetosome-containing digestive vacuoles. Phylogenetic analyses of 16S and 18S ribosomal RNA genes from magnetotactic consortia in marine sediments across the Northern and Southern hemispheres indicate a host-ectosymbiont specificity and co-evolution. This suggests a historical acquisition of magnetoreception by a euglenozoan ancestor from Deltaproteobacteria followed by subsequent diversification. It also supports the cosmopolitan nature of this type of symbiosis in marine anoxic sediments.

MicroScope-an integrated resource for community expertise of gene functions and comparative analysis of microbial genomic and metabolic data.

The overwhelming list of new bacterial genomes becoming available on a daily basis makes accurate genome annotation an essential step that ultimately determines the relevance of thousands of genomes stored in public databanks. The MicroScope platform (http://www.genoscope.cns.fr/agc/microscope) is an integrative resource that supports systematic and efficient revision of microbial genome annotation, data management and comparative analysis. Starting from the results of our syntactic, functional and relational annotation pipelines, MicroScope provides an integrated environment for the expert annotation and comparative analysis of prokaryotic genomes. It combines tools and graphical interfaces to analyze genomes and to perform the manual curation of gene function in a comparative genomics and metabolic context. In this article, we describe the free-of-charge MicroScope services for the annotation and analysis of microbial (meta)genomes, transcriptomic and re-sequencing data. Then, the functionalities of the platform are presented in a way providing practical guidance and help to the nonspecialists in bioinformatics. Newly integrated analysis tools (i.e. prediction of virulence and resistance genes in bacterial genomes) and original method recently developed (the pan-genome graph representation) are also described. Integrated environments such as MicroScope clearly contribute, through the user community, to help maintaining accurate resources.

Evolutionary Evidence of Algal Polysaccharide Degradation Acquisition by Pseudoalteromonas carrageenovora 9T to Adapt to Macroalgal Niches.

About half of seaweed biomass is composed of polysaccharides. Most of these complex polymers have a marked polyanionic character. For instance, the red algal cell wall is mainly composed of sulfated galactans, agars and carrageenans, while brown algae contain alginate and fucose-containing sulfated polysaccharides (FCSP) as cell wall polysaccharides. Some marine heterotrophic bacteria have developed abilities to grow on such macroalgal polysaccharides. This is the case of Pseudoalteromonas carrageenovora 9T (ATCC 43555T), a marine gammaproteobacterium isolated in 1955 and which was an early model organism for studying carrageenan catabolism. We present here the genomic analysis of P. carrageenovora. Its genome is composed of two chromosomes and of a large plasmid encompassing 109 protein-coding genes. P. carrageenovora possesses a diverse repertoire of carbohydrate-active enzymes (CAZymes), notably specific for the degradation of macroalgal polysaccharides (laminarin, alginate, FCSP, carrageenans). We confirm these predicted capacities by screening the growth of P. carrageenovora with a large collection of carbohydrates. Most of these CAZyme genes constitute clusters located either in the large chromosome or in the small one. Unexpectedly, all the carrageenan catabolism-related genes are found in the plasmid, suggesting that P. carrageenovora acquired its hallmark capacity for carrageenan degradation by horizontal gene transfer (HGT). Whereas P. carrageenovora is able to use lambda-carrageenan as a sole carbon source, genomic and physiological analyses demonstrate that its catabolic pathway for kappa- and iota-carrageenan is incomplete. This is due to the absence of the recently discovered 3,6-anhydro-D-galactosidase genes (GH127 and GH129 families). A genomic comparison with 52 Pseudoalteromonas strains confirms that carrageenan catabolism has been recently acquired only in a few species. Even though the loci for cellulose biosynthesis and alginate utilization are located on the chromosomes, they were also horizontally acquired. However, these HGTs occurred earlier in the evolution of the Pseudoalteromonas genus, the cellulose- and alginate-related loci being essentially present in one large, late-diverging clade (LDC). Altogether, the capacities to degrade cell wall polysaccharides from macroalgae are not ancestral in the Pseudoalteromonas genus. Such catabolism in P. carrageenovora resulted from a succession of HGTs, likely allowing an adaptation to the life on the macroalgal surface.

Transcriptome Profiles of Nod Factor-independent Symbiosis in the Tropical Legume Aeschynomene evenia.

Nod factors (NF) were assumed to be indispensable for the establishment of a rhizobium-legume symbiosis until the discovery that certain Bradyrhizobium strains interacting with certain Aeschynomene species lack the canonical nodABC genes required for their synthesis. So far, the molecular dialogue between Aeschynomene and its symbionts remains an open question. Here we report a time course transcriptional analysis of Aeschynomene evenia in response to inoculation with Bradyrhizobium ORS278. The NF-independent symbiotic process was monitored at five time points between bacterial infection and nodule maturity. The five time points correspond to three specific events, root infection by crack entry, nodule organogenesis, and the establishment of the nitrogen fixing process. During the third stage, about 80 NCR-like genes and eight symbiotic genes known to be involved in signaling, bacterial infection or nodulation regulation were highly expressed. Comparative gene expression analyses at the five time points also enabled the selection of genes with an expression profile that makes them promising markers to monitor early plant responses to bacteria. Such markers could be used in bioassays to identify the nature of the bacterial signal(s). Our data represent valuable resources for investigation of this Nod factor-independent symbiosis.

GROOLS: reactive graph reasoning for genome annotation through biological processes.

BACKGROUND: High quality functional annotation is essential for understanding the phenotypic consequences encoded in a genome. Despite improvements in bioinformatics methods, millions of sequences in databanks are not assigned reliable functions. The curation of protein functions in the context of biological processes is a way to evaluate and improve their annotation. RESULTS: We developed an expert system using paraconsistent logic, named GROOLS (Genomic Rule Object-Oriented Logic System), that evaluates the completeness and the consistency of predicted functions through biological processes like metabolic pathways. Using a generic and hierarchical representation of knowledge, biological processes are modeled in a graph from which observations (i.e. predictions and expectations) are propagated by rules. At the end of the reasoning, conclusions are assigned to biological process components and highlight uncertainties and inconsistencies. Results on 14 microbial organisms are presented. CONCLUSIONS: GROOLS software is designed to evaluate the overall accuracy of functional unit and pathway predictions according to organism experimental data like growth phenotypes. It assists biocurators in the functional annotation of proteins by focusing on missing or contradictory observations.